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We extend the hadronic SU(3) non- linear sigma model to include quark degrees of freedom. The 
choice of potential for the deconfinement order parameter as a function of temperature and chemical 
potential allows us to construct a realistic phase diagram from the analysis of the order parameters 
of the system. These parameters are the chiral condensate, for the chiral symmetry restoration, 
and the scalar field $ (as an effective field related to the Polyakov loop) for the deconfinement to 
quark matter. Besides reproducing lattice QCD results, for zero and low chemical potential, we are 
in agreement with neutron star observations for zero temperature. 
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The models used to describe neutron stars can gener- 
ally be divided into two classes. The first class includes 
approaches in which the constituent particles are hadrons 
[lH3| . Some of them incorporate certain symmetries from 
QCD, like chiral symmetry, but they do not include de- 
confinement. Examples of these are hadronic sigma mod- 
els (343- The second class includes quark star models, 
which usually do not directly incorporate hadronic de- 
grees of freedom in the model formulation. Examples 
of these are bag-model studies ja as well as quark-NJL 
model and quark sigma- models |9[. 

Using these approaches hybrid neutron stars, which 
consist of a hadronic and a quark phase, are normally 
described by adopting two different models with sepa- 
rate equations of state for hadronic and quark matter 
(see e.g. [iol|). They are connected at the chemical po- 
tential in which the pressure of the quark EOS exceeds 
the hadronic one, signalling the phase transition to quark 
matter. Within our approach we employ a single model 
for the hadronic and for the quark phase. 

The extension of the hadronic SU(3) non-linear sigma 
model to quark degrees of freedom is constructed in a 
spirit similar to the PN JL model [ll| , in the sense that it 
is a non-linear sigma model that introduces a scalar field 
which suppress the quark degrees of freedom at low den- 
sities/temperatures. In QCD this scalar field was named 
Polyakov loop and is defined via $ = |Tr[exp (i J drA^j], 
where A4 — iAo is the temporal component of the SU(3) 
gauge field. In our case, this scalar field is also called $, 
in analogy to the PNJL approach with an effective po- 
tential for the field, as discussed below, that drives the 
phase transition in the field $ representing a phenomeno- 
logical description of the transition from the confined to 
the deconfined phase. 

The Lagrangian density of the non-linear sigma model 
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in mean field approximation reads: 

L = LjCin + Lint + Lself + LsB — U, (1) 

where besides the kinetic energy term for hadrons, 
quarks, and leptons (included to insure charge neutrality) 
the terms: 



Lint 



(2) 



-Self 
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L SB = mlUcr + (V2m 2 f k - -±=mlU) C, 
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represent the interactions between baryons (and quarks) 
and vector and scalar mesons, the self interactions of 
scalar and vector mesons and an explicit chiral symme- 
try breaking term, responsible for producing the masses 
of the pseudo-scalar mesons. The $ potential U will be 
discussed in the following. The underlying flavor sym- 
metry of the model is SU(3) and the index i denotes the 
baryon octet and the three light quarks. The mesons 
included are the vector-isoscalars u> and <f>, the vector- 
isovector p, the scalar-isoscalars a and C (strange quark- 
antiquark state) and the scalar-isovector S. The isovector 
mesons affect isospin-asymmetric matter and are conse- 
quently important for neutron star physics. The coupling 
constants of the model are shown in Table. [TJ They were 
fitted to reproduce the vacuum masses of the baryons and 
mesons, nuclear saturation properties (density po = 0.15 
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TABLE I: Coupling constants for the model containing only 
baryons 



g Nw = 11.90 

QNcr = -9.83 

<? Aw = 7.93 

QAcr = -5.52 

fc = 1.19 
fc 3 = 2.65 



gN<p = 
g NS = -2.34 
gA4> = -7.32 

gAs = 
fci = -1.40 
k 4 = -0.06 



g Np = 4.03 
g N( = 1-22 

9a p = 
.9 Ac = -2.30 
k 2 = 5.55 
g 4 = 38.9 
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FIG. 1: (Color online) Order parameters for chiral symmetry 
restoration and deconfinement to quark matter as a function 
of baryonic chemical potential for star matter at zero temper- 
ature. 



fm 3 , binding energy per nucleon B/A = —16.00 MeV, 
nucleoli effective mass = 0,67 Mat, compressibility 
K = 297.32 MeV), asymmetry energy (E sym = 32.50 
MeV), and reasonable values for the hyperon potentials 
(U A = -28.00 MeV, Us = 5.35 MeV, U s = -18.36 MeV). 
The vacuum expectation values of the scalar mesons are 
constrained by reproducing the pion and kaon decay con- 
stants. A detailed discussion of the purely hadronic part 
of the Lagrangian can be found in 0, [H, HI]- 

The mesons are treated as classical fields within the 
mean- field approximation [l4|. Finite-temperature cal- 
culations include the heat bath of hadronic and quark 
quasiparticles within the grand canonical potential of the 
system. It is the defined as: 



— —Lj nt 

k r 



LSelf — L SB — Ly ac 

^E^f - d^Miie-^rw-tf)), (5) 

where Ly ac in the vacuum energy, ji the fermionic degen- 
eracy, E*(k) = \Jk 2 + M* 2 the single particle effective 
energy, and n\ = /i; - g iu] u> - g^ - g ip r 3 p the effective 
chemical potential of each species. The chemical poten- 
tial for each species //$ comes from the chemical equilib- 



rium conditions. Finite temperature calculations also in- 
clude a gas of free pions and kaons. As they have very low 
mass, they dominate the low density/ high temperature 
regime. All calculations were performed considering zero 
net strangeness except the zero temperature star mat- 
ter case since, for neutron stars, the time scale is large 
enough for strangeness not to be conserved. 

The effective masses of the baryons and quarks are 
generated by the scalar mesons except for a small explicit 
mass term Mq (equal to 150 MeV for nucleons, 354 MeV 
for hyperons, 5 MeV for up and down quarks and 150 
MeV for strange quarks) and the term containing $: 

M* B = g Ba a + g BS T 3 S + g B (( + M 0b + g B ^<i> 2 , (6) 



M q = 9qoO + g q sT 3 S + g qC ( + M Qq + g q9 (l - $). 



(7) 



With the increase of temperature/density, the a field 
(non-strange chiral condensate) decreases its value, caus- 
ing the effective masses of the particles to decrease to- 
wards chiral symmetry restoration. The field $ assumes 
non-zero values with the increase of temperature/density 
and, due to its presence in the baryons effective mass 
(Eq. |B])), suppresses their presence. On the other hand, 
the presence of the $ field in the effective mass of the 
quarks, included with a negative sign (Eq. ([7])). in- 
sures that they will not be present at low tempera- 
tures/densities. As can be seen from the different orders 
of $ and different signs in the new effective mass terms, 
the motivation for this construction is not derived from 
QCD. It is a simple effective way to change degrees of 
freedom within the same model. Note that in the PNJL 
approach the coupling of the quarks to the Polyakov loop 
can be derived to be included in the quark and antiquark 
distribution functions in the grand canonical potential. 
However, this leads to non-vanishing quasi-quark contri- 
butions at any temperature below T c , which we avoid in 
our phenomenological approach (Eqs. (6,7)). 

The behavior of the order parameters of the model is 
shown in Fig. [1] for neutron star matter at zero tempera- 
ture. The difference between this kind of matter and the 
so-called symmetric matter comes from the assumption 
of charge neutrality, essential for the stability of neu- 
tron stars, and beta equilibrium. In this case, the chi- 
ral symmetry restoration, which is a crossover for purely 
hadronic matter, turns into a first order phase transi- 
tion by the influence of the strong first order transition 
to deconfined matter. The model is consistent in the 
sense that both order parameters are related. The small 
increase in the chiral condensate value during the tran- 
sition is due to the smaller quark baryon number (1/3) 
compared to the baryonic one. 

The effective normalized masses of baryons and quarks 
show the relation between this quantities and the order 
parameters, responsible for the dynamics of the model 
(Fig. [T]and[2]). Since the coupling constants in the $ 
term of the effective mass formulas are high but still fi- 
nite, the effective masses of the degrees of freedom not 
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FIG. 2: (Color online) Effective normalized mass of different 
species as a function of baryonic chemical potential for star 
matter at zero temperature. 
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FIG. 3: (Color online) Phase diagram: temperature versus 
baryonic chemical potential. The lines represent first order 
transitions. The circles mark the critical end-points. 



TABLE II: Additional coupling constants for the model con- 
taining baryons and quarks 



9qu = 

g qa = -3.00 
a = -1.85 
a 3 = -0.40 



9q4> = 
9qS = 

o a = -1.44xl0~ 3 
gs* = 1500 MeV 



9qp = 

g qC = -3.00 
a 2 = -0.08 
3,2* = 500 MeV 



To = 200 MeV 



To = 270 MeV for pure gauge case 



effectively present in each phases are high but also finite. 
The effective masses normalized by the baryonic number 
are shown in Fig. [21 These quantities are directly related 
to the onset of particles appearance in the system, that 
reads: 



M* 



^ e Qb. 



(8) 



where Qg { is the baryonic number, fi e is the electron 
chemical potential, Qi is the electric charge, lis the 
strange chemical potential, and Qsi the strangeness of 
each species. 

Continuing the analogy to the PNJL model, the po- 
tential U for $ reads: 



U = (a T 4 + a lL L A 
+a 3 T 4 log(l -6$ 2 



8$ 3 -3$ 4 ). 



(9) 



It is a simplified version of the potential used in [TH, [T(| 
and adapted to also include terms that depend on the 
chemical potential. The two extra terms (that depend 
on the chemical potential) are not unique, but the most 
simple natural choice whose parameters are chosen to re- 
produce the main features of the phase diagram at finite 



densities. While in the first part of the potential $ 2 en- 
sures U to be bound from below, in the second part of 
the potential T 4 log (1 - 6$ 2 + 8$ 3 - 3$ 4 ) ensures that 
$ is always (for any region of the phase diagram) bound 
between zero and one. 

The coupling constants for the quarks are shown in 
Tab. [IT] and are chosen to reproduce lattice data as well 
as known information about the phase diagram. The 
lattice data includes a first order phase transition at 
T = 270 MeV and a pressure functional P(T) similar 
to Refs. [H, [l(| at Li = for pure gauge, a crossover at 
vanishing chemical potential with a transition tempera- 
ture of 171 MeV (determined as the peak of the change 
of the chiral condensate and $) and the location of the 
critical end-point (at li c = 354 MeV , T c = 167 MeV for 
symmetric matter in accordance with one of the existent 
calculations [l7|). The phase diagram information in- 
cludes a continuous first order phase transition line that 
terminates on the zero temperature axis at four times 
saturation density. 

As can be seen in Fig. [31 the transition from hadronic 
to quark matter obtained is a crossover for small chemical 
potentials. Beyond the critical end-points a first order 
transition lines for symmetric as well as for star mat- 
ter begin. The critical temperatures for chiral symmetry 
restoration coincide with the ones for deconfmement in 
both cases. Since the model is able to reproduce nuclear 
matter saturation at realistic values for the saturation 
density, nuclear binding energy, as well as the compress- 
ibility and asymmetry energy, we also show a line in the 
phase diagram for the nuclear matter liquid-gas phase 
transition. 

One way to test the model and to compare its pre- 
dictions with known observational data is to study the 
high density/low temperature part of the phase diagram 
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FIG. 4: (Color online) Population (baryonic density for dif- 
ferent species as a function of baryonic chemical potential) for 
star matter at zero temperature using local charge neutrality 
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FIG. 5: (Color online) Mass-radius diagram 



and compare our results with neutron star observations. 
The critical point for star matter lies at a slightly higher 
chemical potential than for the symmetric case and the 
first order transition line terminates on the zero temper- 
ature axis at fis — 1345 MeV (grey line in Fig. Up 
to this point, the charge neutrality was considered to be 
local, meaning that each phase had to be charge neutral 
by itself. At finite temperature the two phases contain 
mixtures of hadrons and quarks, which are dominated by 
hadrons or quarks, depending on the respective phase. At 
vanishing temperature there is no mixture, i.e. the sys- 
tem exhibits a purely hadronic and purely quark phase 
(Fig. 5]) . The density of electrons and muons is significant 
in the hadronic phase but not in the quark phase. The 
reason for this behavior is that because the down and 
strange quarks are also negatively charged, there is no 
necessity for the presence of electrons to generate charge 
neutrality, and only a small amount of leptons remains 
to assure beta equilibrium. 

The quarks are totally suppressed in the hadronic 
phase and the hadrons are suppressed in the quark phase 
until a certain chemical potential (above 1700 MeV for 
T = Q). This behavior comes from the fact that the 
coupling constants in the $ term of the effective mass 
formulas are high but still finite, so at very high chem- 
ical potential the threshold in Eq. ([SJ can be reached a 
second time for hadrons. This threshold, that is higher 
than the density in the center of neutron stars, estab- 
lishes a limit for the applicability of the model. The hy- 
perons, in spite of being included in the calculation, are 
suppressed by the appearance of the quark phase. Only 
a very small amount of A appears right before the phase 
transition (Fig. 0]) . The strange quarks appear after the 
other quarks and also do not make substantial changes 
in the system. 

The possible neutron star masses and radii are cal- 
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FIG. 6: (Color online) Population (baryonic density for dif- 
ferent species as a function of baryonic chemical potential) for 
star matter at zero temperature using global charge neutral- 
ity. 



culated solving the Tolmann-Oppenheimer-Volkof equa- 
tions [ TsL figj l . The solutions for hadronic (same model 
but without quarks) and hybrid stars are shown in Fig. [SJ 
where besides our equation of state for the core, a sep- 
arate equation of state was used for the crust [20]. The 
maximum mass supported against gravity in our model 
is 2.1M Q in the first case and around 2.0M Q in the sec- 
ond. Because the equation of state for quark matter is 
much softer than the one for hadronic matter, the star 
becomes unstable right when the central density is higher 
than the phase transition threshold. 

There is still another possible option for the configu- 
ration of the particles in the neutron star [21]. If instead 
of local we consider global charge neutrality, we find a 



5 



— local charge neutrality 




e (MeV/fm 3 ) 



FIG. 7: (Color online) Equation of State (pressure as a func- 
tion of energy density) for star matter at zero temperature 
using local and global charge neutrality. 

mixture of phases. This possibility, which is a more re- 
alistic approach, changes the particle densities in the co- 
existence region making them appear and vanish in a 
smoother way (Fig. [5]). Therefore, the maximum mass 
allowed for the star is slightly lower in this case than in 
the previous one, as can be seen from the dotted line 
in Fig. however, this possibility allows stable hybrid 
stars with a small amount of quarks. The mixed phase 



constitutes the inner core of the star up to a radius of ap- 
proximately 2km. The equation of state for both cases is 
shown in Fig. [7] The large jump in the pressure in the lo- 
cal charge neutrality case explains why the neutron stars 
become immediately unstable after the phase transition 
in this configuration. 

We conclude that our model is suitable for the descrip- 
tion of neutron stars. The maximum mass predicted is 
around the most massive pulsars observed [22-26]. The 
radii lie in the allowed range being practically the same 
for hadronic or hybrid stars. A major advantage of our 
work compared to other studies of hybrid stars is that 
because we have only one equation of state for differ- 
ent degrees of freedom we can study in detail the way 
in which chiral symmetry is restored and the way decon- 
fmement occurs at high temperature/density. Since the 
properties of the physical system, e.g. the density of par- 
ticles in each phase, are directly connected to the order 
parameter for deconfinement $ it is not surprising that 
we obtain different results in a combined description of 
the degrees of freedom compared to a simple matching of 
two separate equations of state. 

Since the model additionally shows a realistic structure 
of the phase transition over the whole range of chemical 
potentials and temperatures as well as phenomenologi- 
cally acceptable results for saturated nuclear matter, this 
approach presents an ideal tool for the study of ultrarel- 
ativistic heavy-ion collisions. Calculations along this line 
are in progress (27| . 
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